Universality of quantum critical dynamics in a planar OPO 
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We analyze the critical quantum fluctuations in a coherently driven planar optical parametric 
oscillator. We show that the presence of transverse modes combined with quantum fluctuations 
changes the behavior of the 'quantum image' critical point. This zero-temperature non-equilibrium 
quantum system has the same universality class as a finite-temperature magnetic Lifshitz transition. 
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The non-equilibrium system consisting of a nonlinear 
crystal inside a laser driven Fabry-Perot interferometer, 
that couples sub-harmonic intra-cavity modes to a har- 
monic pump pj, is known as an optical parametric oscil- 
lator (OPO). As well as demonstrating quantum squeez- 
ing |3] and EPR entanglement in numerous quantum in- 
formation experiments, this system is widely used in fre- 
quency conversion applications. In single- mode experi- 
ments, there is a critical point in the phase diagram. This 
is caused by an increase in the pump intensity, which 
results in a transition from a disordered (but quantum 
squeezed) phase below threshold, to an ordered phase 
with a coherent output above threshold. 

When extended to an interferometer with multiple 
transverse modes, more complex dynamical effects oc- 
cur due to diffraction of the down-converted light, 
which are governed by the Swift-Hohenberg equation 
near threshold^. The theory can be quantized^, and 
hence includes quantum fluctuations^]. Experimentally 
observable@,0] 'quantum images' are evidence for quan- 
tum pattern formation with spatio-temporal correlations 
in the output quadratures. Thus, the system can show 
both spatial critical fluctuations and non-equilibrium 
spontaneous pattern formation, which occurs in many 
fields of physics and other sciences ja]. 

In this Letter we apply the theory of finite size scaling 
to solve for the critical quantum dynamical properties, 
and obtain the universality class this phase-transition 
corresponds to. We provide a full quantum descrip- 
tion of this non-equilibrium system using the positive 
P-representation, focusing on the nature of the critical 
point and critical fluctuations. This allows us to obtain 
an analytic solution for the functional distribution of the 
large critical fluctuations caused by quantum noise in the 
down-conversion process. There is an unexpected uni- 
versality property in the solutions. Even though this is a 
non-equilibrium quantum system of coupled boson fields, 
we find the in one quadrature, the quantum fluctuations 
have exactly the same behavior as a classical thermal sys- 
tem of fields at a two-dimensional Lifshitz point, which is 
a model commonly used to describe the phase transition 
to a modulated magnetic phase; in the complementary 
quadrature there is strong entanglement. 



The unitary evolution of the OPO system can be de- 
scribed by the Hamiltonian0] 



* = £ 



n=0,l 



d 2 x I At 



2uj„ 



A n 



-ill I d 2 xA \£*e 2lULt 



X 



A\ 2 } 



(1) 



The term x is a coupling parameter that depends on 
the nonlinear crystal, the frequencies of the field modes 
are wi, u>q = 2u>i, v is the intracavity group velocity, and 
A n is the n — th photon field. The pump is described by 
the amplitude £ that could carry a spatial structure - but 
here we will assume a constant plane wave input. In ad- 
dition, there are damping effects due to output couplings 
from the cavity mirrors, which can be well approximated 
using as a Markovian master equation for the density 
matrix p, so that: 
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where £„ [p] = J d 2 x 2A n pA n - pA n A n - A n A n p 
scribes the output coupling from the nth intra-cavity 
mode, with damping rate 7„. This leads to a set of 
Fokker-Planck equations, mapped from the the quan- 
tum density matrix, using operator representation the- 
ory. These are valid provided boundary terms vanish in 
the mapping transformation, which we have checked nu- 
merically. Using the positive P-representation, we derive 
the following stochastic equations^, 0, 0] in a rotating 
frame at frequency ujl- 
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Here we write the two dimensional Laplacian causing 
diffraction, as V 2 = d 2 /dx 2 + d 2 /dy 2 . The complex 
relaxation rates are ji = 7,(1 + iAi). The relative de- 
tunings between the pump laser at 2ujl, and the modes 
supported by the cavity are Aq = (ujq — 2ujl)/jo, and 
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Ai = (lui — ujl)Hi • The diffraction rate is defined as 
D = w 2 /(27iwi). The stochastic field £1 which describes 
quantum noise is real and Gaussian, with correlations of 
(&(t)> =0and (£i(a?,t)£i(a? ) f)) = 5 2 {x - x^)8{t - t'). 

In addition, there are equations that correspond to 
the hermitian conjugate fields. As elsewhere in this Let- 
ter, we obtain these by conjugating the constant terms 
and replacing the stochastic and noise fields accord- 
ing to: Ai — > Af, £1 — > £y, where £1, are inde- 
pendent real Gaussian noises. These two noise fields 
are sufficient to generate all quantum effects, and are 
physically caused by the discrete nature of the pho- 
ton pairs produced in down-conversion. The c-number 
fields Ai(t,x), Af(t,x) are therefore not complex con- 
jugate, although they are stochastically equivalent in 
terms of normally-ordered operator moments to pho- 
ton operator fields Ai(t,x), A\{t,x). Thus for exam- 
ple, the photon number density is(Jlt (t, x)Ai(t, a?)) = 
(A+(t,x)Mt,x<)} . 

If we remove the transverse modes from the above 
equations we return to the well known single mode OPO 
theory. This system has a quantum critical point - a 
phase transition in the infinite volume limit, where the 
quantum fluctuations are reduced below the vacuum level 
for the squeezed quadrature and become huge for the un- 
squeezed quadrature 0. In a recent analysis 0] of this 
problem near the critical point, going beyond the lin- 
ear theory, we obtained a scaling law for the squeezing 
quadrature spectrum near threshold, and the parameters 
for the optimum squeezing. 

The introduction of transverse modes generates a spa- 
tial structure in the sub-harmonic field with an inten- 
sity correlation function known as a "quantum image" 0, 
since it is supported by quantum fluctuations. To treat 
this problem analytically, we can perform an adiabatic 
elimination of the stable pump mode in the limit of 
7o 7i and A — > 0. That is, we assume that the 
pump mode has a short relaxation time. 

Neglecting pump diffraction - which is negligible in the 
critical regime - we obtain an adiabatic solution for the 
pumped field: Aq = {£ — x*A\/i) /70, together with a 
similar equation for the conjugate term. This solution 
takes into account the depletion of the pumping mode 
that supplies energy for down-converted light, and leads 
to an adiabatic equation for the down-converted field: 



dA x 
dt 



(4) 



Next, we introduce dimensionless variables r = t/to 
and f= x/xq, with a corresponding down-converted field 
a = xqAi. Due to critical slowing down, the character- 
istic length xq and time to scale as to = l/(<7 7 i) and 
Xq = D j J~g~ c , where the effective nonlinear coefficient is 



9c = |x| 4 ^ 3 /[8£ ) 7o7i] 2 ^ 3 , and we will assume that j c « 1. 
The dimensionless driving field is Jl = X^/[7i7o] = H+iO- 
We also introduce appropriately scaled noise fields with 
£(r,r) = x \/tocli(t,x), and the corresponding hermitian 
conjugate terms. 

With these definitions, we find that: 
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This equation includes both the down-conversion term 
proportional to , which generates a squeezed signal — 
together with a nonlinear saturation term proportional 
to — a\a\ , which limits the down-converted amplitude, 
and leads to finite size critical fluctuations. 

Critical fluctuations are most usefully analyzed with 
scaled quadratures that correspond to experimentally ac- 
cessible homodyne detection. These are defined as 



Y(t,t) = i \ciy (r, r) — ot\ (r, r)l . 
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Similarly, there are quadrature noise fields defined 
as £r(T,r) = f(r,r) + £+(V,f), and t] v (t,F) = 
i\/9c [£ + ( T : r) — £(t, r)]. The resulting quantum dynam- 
ical equations for these signal field quadratures are: 
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The linear decay matrix that couples the X and Y 
quadratures is given by: 
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We assume also that close to threshold, 
small enough detunings, j x — 0(1) and j xy = — (0 + 
Ai)/y^fc = 0(1). We can always choose quadrature 
phases so that 9 = Ai + 0(g c ). With this choice, Y is 
mainly coupled to the X quadrature via the diffraction 
term, which couples noise from the critical fluctuations 
back into the squeezed quadrature. This implies that 
jl = 1 + 0(g c ) and 7y = 2 + 0(^/g2), so that the noise 
correlations are given by; 

(Ur^Ur'^)) = 2^(r - r')5 2 (r - ?) + 0(g c ) . (9) 

We can now perform a second type of adiabatic elimi- 
nation, which is valid in a neighbourhood of the critical 
point. This takes into account the fact that the fluctua- 
tions in the X quadrature become very slow near thresh- 
old, while the Y quadrature still responds on fast time- 
scales of order I/71 . To leading order we can drop terms 
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of 0(y/g2) where g c <C 1, and approximate the above 
equations as follows: 



^ = - lx X- lxy Y-X 3 -V 2 Y- 







-2Y + V 2 A . 
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We can therefore eliminate the fast or non-critical 
quadrature variable Y, by writing the steady state solu- 
tion of the Y quadrature as Y ~ V 2 A/2 . This produces 
a reduced equation for the critical quadrature variable 
X, which is valid near threshold: 



dX 
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The above Langevin equation is a Ginzburg-Landau 
equation describing the critical quadrature dynamics. 
Unlike the usual application of this equation, we note 
that system is a non-equilibrium one. The noise term 
£ x is of quantum origin rather than thermal origin, and 
is present at zero temperature. It is possible to write 
an equivalent functional Fokker Planck equation for the 
probability density P[-X], 
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and look for the equilibrium distribution in the form 
P[X] = Nexp(-V[X]), where V(X) is a potential func- 
tional. Making this substitution, the solution for the dis- 
tribution P[X] is given by : 
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This expression is exactly the same as the Ginzburg- 
Landau free energy of a next nearest neighbor magnetic 
interaction, where X plays the role of an order parame- 
ter. That is, we have been able to map this problem into 
a soluble magnetic phase-transition equation with a Lif- 
shitz point [11]. The phase diagram of this optical system 
should therefore have two ordered phases, one of them a 
spatially modulated phase associated with a pattern for- 
mation. This generic behavior is known to occur in an 
OPO, from previous analysis jj. 

In this analogy, the "optical paramagnetic phase" cor- 
responds to a random photon emission from the OPO 
operating below threshold, the "optical ferromagnetic 
phase" to a continuum emission uniformly distributed 
in the transverse plane parallel to the cavity operating 
above threshold. In the "optical ferromagnetic modulated 
phase", we have a continuum emission but with modu- 
lated quadrature in this plane. At the Lifshitz point, all 
three phases co-exist. 

The line 7^ = is the line of the second-order phase 
transition between order-disorder (coherent-incoherent) 
states. In the incoherent phase below threshold, j x > 0, 
and in the uniform coherent phase above threshold j x < 




Figure 1: Variance of the quadrature Fourier component 
(11) i\X(k)\ 2 ) as a function of y x . Results obtained on a 100 x 100 
lattice with a 40 x 40 domain size, using periodic boundary 
conditions and averaging over 100 stochastic trajectories. 



0, as expected in the single mode case. If ~f x vanishes 
we have a Lifshitz point over the line j x = 0, and thus 
a triple point characterizing the coexistence of the three 
phases. This holds in the case of perfect tuning of the 
signal field inside the cavity, so that Ai = 0. 

In condensed matter physics, the nature of the Lif- 
shitz point 0, Q] is crucially dependent on the order 
parameter and spatial dimension. Depending on the di- 
mensionality of the order parameter, the system may or 
may not have a true phase transition in two dimensions. 
Our system has a one-dimensional real order parame- 
ter and two spatial dimensions with transverse modes, 
so we expect that this system should have a true phase 
transition in the infinite volume limit at finite temper- 
ature. According to the Mermin- Wagner theorem yj| a 
increasing the order parameter dimension (as in type II 
down-conversion) would result in phase fluctuations that 
completely destroy any long range order. Similarly, re- 
ducing the spatial dimension would result in a continuous 
transition without a threshold. 

A numerical simulation of the Ginzburg-Landau equa- 
tions ifTT)) with a continuously scanned input shows that 
the sub-harmonic quadrature correlations appear to have 
a true critical point in two transverse dimensions. This 
result is shown in Figure (1), which graphs (|A(fc)| 2 ) 
around k — 0, as a function of the driving field j x near 
threshold. 

While these large fluctuations are occurring, we note 
that there are still strong non-classical correlations in the 
squeezed quadrature. This can be seen by analysing the 
relevant equations to the next order in g c , which we also 
simplify by using a Gaussian factorization [jg : 



BY 

9c-^ « -lyY + V 2 A - g c Y{X 2 ) - % X X + { y , (14) 

where ^ yx = j yx + 2g c (XY). In general, one can always 
choose an optimum local oscillator phase so that 6 = 
k 2 + Ai + 2g c (XY), in order to minimise the feedback 




Figure 2: Squeezing spectrum as a function of frequency, with 
(solid line) and without (dotted line) nonlinear corrections. 
The parameters used are 7^ = 0.5, g c = 0.01. 



of critical fluctuations into the squeezed quadrature at 
a given transverse momentum k. This leads to Fourier 
solutions which showing that entanglement between 
the modes of momentum k and — k can still occur at small 
enough wave- vectors, resulting in a universal squeezing 
spectrum as a function of frequency: 



v(n) = 1 



1 - g c ((X 2 ) + lx ) 



(g c n/2y + (l+g c ((X^- lx )/2Y 



(15) 



This result differs from the linearised predictions of 
earlier treatments 0. A graph of the resulting spectrum 
in the Gaussian approximation is shown in Fig (2) , com- 
pared to the linearized squeezing spectrum, showing large 
differences near threshold. 

In summary, we have shown that the planar non- 
equilibrium OPO with quantum noise can be mapped to 
a magnetic phase-transition in two dimensions. Since the 
present case has a scalar, real order parameter it is anal- 
ogous to the uni-axial (to = 1) magnetic order parame- 
ter case, which is known to have a thermal equilibrium 
Lifshitz-point phase-transition at finite temperature^]- 
This demonstrates a striking resemblance between known 
thermal equilibrium phase-transitions, and a quantum 
non-equilibrium system in which quantum noise replaces 
thermal noise. In this system there are also quantum 
correlations of the emitted photons, causing quantum 
squeezing and entanglement. Nevertheless, this highly 
non-classical behavior is found only in the squeezed (Y) 
quadrature which has no critical slowing down - and co- 
exists with a rather classical and universal critical fluc- 



tuation field in the conjugate (X) quadrature. 
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